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MEMORANDUM 


TO: 


N.  Peacock  and  the  Culham  Group 
(Memorandum  Number  Two) 


FROM:  M.  V.  Goldman 


SUBJECT:  Linear  Inhomogeneous  Theory  of  CO^Enhanced  Langmuir 


Turbulence 


I.  The  scattering  cross  section  for  the  plasma  line  is  given  roughly  by 
2 


O = -T  N 6!)  7 , 
2 e s 
a 


whe re. 


/ 2 \ 
\m  c / 


•,  „ ,„-26  2 

7.9  x 10  cm  , 


2 V 2/4wi2  ^sX"1 

“ =~2  ) 
k ' c ' 


6 = fid)  cos(0  - 60  ) - 6<J>  cos(0  + 60  ) , 
s s s s s 


J = lim 


<|6E(ky)|  > 


v*®°  4TT0  V 
e 


(r-?S  u. 

\ J VT  ~*’ao  41T0  VT  / 


(1) 


(2) 


(3) 


(4) 


(5) 


2 


In  this  expression,  N is  the  total  number  of  scattering  electrons,  which 

e 

is  equal  to 

N = n ir  2$  6r  , (6) 

e o o c 

19  -3 

where  n^  e 10  cm  is  the  density  of  electrons  in  the  neighborhood 

where  a)  = u)  , 1 is  the  axial  size  of  the  ruby  spot  (smaller  than  the 
C02  pe 

C0„  spot),  r is  the  radius  at  which  n(r  ) ~ n , <t>  is  the  maximum  azimu- 
2 o o o c 

thal  angle  for  which  the  CO^  beam  reaches  the  frequency  matched  radius: 


and  (5r  is  the  radial  thicJcness  over  which  I is  close  to  its  maximum. 


II.  For  sinplicity  in  treating  plasma  wave  convection  effects,  we  shall 
replace  the  volume  in  Fig.  1 by  a slab,  as  illustrated  in  Figure  2: 


3 


N I 


Vn 


^C0„ 


Note  that,  in  this  geometry  the  Langmuir  waves  probed  by  the  ruby  laser 
(driven  by  ) are  propagating  in  the  y-direction.  Inhomogeneity  is 

allowed  only  in  the  x-direction. 

III.  The  spectrum,  I , of  Langmuir  waves  excited  by  the  electron-ion 
decay  instability  is  the  solution  to  a WKB  kinetic  equation: 


— + — - — + — - - 2 (Y  - Y ~ YNL({I})  ) I ~ 2S  (8) 

at  3k  3x  3x  3k  (Yg  Yc  Y (t  i,)  I - 2S  . (8) 

x x 


Here,  we  have  assumed  a convective  instability , for  which  the  growth 
rate  is  y . Yc  is  the  linear  collisional-damping  rate.  S is  an  appro- 
priate  spontaneous  emission  term,  and  y is  a possible  nonlinear  satura- 
tion term  which  would  be  a functional  of  J.  In  memo  1,  we  estimated 
6x  and  I by  ignoring  x and  t-dependence  in  (8) , and  using  a nonlinear ly 
saturated  I.  In  this  memo,  we  will  ignore  the  nonlinear  term  yNL  and 


find  6x  and  I by  finding  a steady  state  inhomogeneous ly  saturated  J. 
This  is  in  the  spirit  of  Perkins  and  Flick,  but  has  physical  effects 


4 


ignored  by  them.  Our  treatment  is  WKB.  If  x = x(t)  and  k = k (t) 

x x 

define  the  WKB  orbits  of  the  Langmuir  waves,  then  along  these  orbits, 

T(k  (t)  , x(t)  ; k , <i»  ) obeys 
x y L 


-TZ  I - 2 (Y  - Y ) I = 2S(k  (t),  x ( t) ) • 
dt  g c x 


In  steady  state. 


d ^ x 3 , 3x  a 
dt  3t  3t  3t  3t 


The  WKB  orbits  are  the  solution  to  Hamilton's  equations 


3k  3u)  _ 3(u 

x L 3x  _ L 

3t  3x  ' 3t  3k 


u.  = (w  2(x)  + 3v  2[k  2 + k 2])1/2 
L p e 1 x y J 


We  assume  a linear  density  profile: 


3u)  (i) 

3x  2L 


(ignore  dependence  on  x on  right-side) 


The  solutions  to  (11)  and  (12)  are  then 


V 

kx  = - it  ' 


, v 2t2  3Lk  2 

3 e x 

x = x - — = X — , 

O 4 L o ,2 
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L 


where  is  the  reflection  point  of  the  Langmuir  wave.  In  addition 

to  (14)  and  (15)  we  have  the  constants  of  the  motion,  k^  and  Thus, 

x is  the  solution  to 
o 


2 2 2 
(0)  (x  ) + 3v  k ) . 
p o e y 


(16) 


We  must  follow  a Langmuir  wave  from  a point  x < x to  x where  it  has  its 

o o 

k in  the  y-direction,  and  can  be  probed  by  the  ruby.  The  weakness  of 
our  theory  is  that  WKB  tends  to  break  down  at  the  reflection  point. 
Nevertheless,  we  have  shown  (see  Appendix)  that  this  same  WKB  method 
can  predict  the  Perkins-Flick  threshold  to  within  a factor  of  two,  under 
the  relevant  conditions  (T  » T J for  their  non-WKB  theory.  Figure  3 
shows  the  trajectory  of  a Langmuir  wave  which  can  be  probed  by  the  ruby: 


FIGURE  3 


We 


wish  to  study  orbits  for  which  the  spectrum,  I,  is  a maximum  at  x = x 


and  is  one-half  of  its  maximum  at  x = x - 6x.  This  will  determine  6x. 

o 

(This  is  a different  criterion  from  memo  1.) 


The  kinetic  equation,  (9)  may  be  integrated  to  give 

r 0 -2)Voat  1 

lit  = 0)  = e I 1 + / dt'  e i S(t’)  I , 


A = / dt  fg(t)dt  , 


r (t)  = y it)  - y 

g g c 


Here,  we  are  letting  t = 0 be  the  time  at  which  reflection  occurs 

(x  = x , k =0),  and  t.  < 0 be  some  initial  time  (see  Fig.  3).  We 
oxi 

have  assumed  I = 1 at  t = t^,  corresponding  to  an  equilibrium  spectrum. 
The  dominant  contribution  to  the  spontaneous  emission  comes  from  the 
"beating"  of  low- frequency  ion  Cerenkov  emission  with  the  CO^  electric 
field,  allowing  it  to  be  up-shifted  into  a source  for  Langmuir  waves.  It 
can  be  shown  (DuBois  and  Goldman,  Phys.  Rev.  164,  207  (1967))  that: 


S = Y ♦ c — 

C 6d  -fci 

o L 


The  first  term  is  Bremsstrahlung-emission  of  Langmuir  waves,  and  the 

second  is  "beat-emission,"  which  is  proportional  to  the  CC>2  pump  intensity 

through  Y • If  we  write  y = T + y > the  dominant  contribution 
g g g c 

from  the  right-side  of  (17)  is  from  the  integral  over  T uj  / (ou  -ojt  ) , 

g p o L 

which  may  be  performed  exactly: 


I it  = 0)  = ^ — e 

U)  “(it 
O L 


(eI 2*  » 1) 


•a 
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Now  it  remains  to  do  the  integral  A given  in  equation  (18).  For  this 


we  require  an  expression  for  the  growth  rate  Y . When  Y <<  w , the 

g g a 


growth  rate  is  given  by 


(k*E  ) u> 

° _E 


32TTnO 


Im  [(k2AD2)  £ (k.  wL~Wo)j 


-1 


(22) 


where  £ is  the  lew  frequency  dielectric  function,  and  is  the  local 


value  of  C0_  electric  field.  When  T >>  T.f  ion-acoustic  waves  are 
2 e i 


long-lived,  and  a resonant  approximation  for  £ may  be  used.  Then  (22) 
and  (18)  may  be  combined  to  give  the  Perkins-Flick  "threshold"  (see 


appendix).  However,  when  T - T. , as  in  this  experiment,  one  cannot  use 

e l 


a resonant  approximation.  In  fact,  the  instability  should  really  be 
thought  of  as  the  decay  of  a pump  photon  into  a Langmuir  wave  plus 
discrete  ions.  Rather  than  use  the  exact  £ (see,  for  example,  DuBois 
and  Goldman,  Phys.  Rev.  Lett.  544  (1965)),  an  adequate  approxima- 

tion at  equal  temperatures  is 


2 2 2 

(D  -a)  - (u)  -u)  ) /2k  c 

r 2 _ 2S  -i-l o L _ o L s 

Im  [jk  Ad  ) £ J 


kc 


(23) 


1/2 


c = 1-7  v.  = 1.7  (m  /m.)  v , 
s l e i e 


(24) 


-1/2 

This  "resonant  function"  has  a maximum  value  of  e =0.61  when 


cu  - w.  (k)  = kc 
O L S 


(25) 


which  therefore  constitutes  the  frequency-matching  condition.  Both 
the  maximum  value  and  (25)  are  very  close  to  the  exact  values. 

We  shall  include  another  effect  missing  from  the  Perkins-Flick  analysis: 
the  swelling  of  the  CO^  electric  field  due  to  geometric  optics.  If 
IVAC  is  the  incident  vacuum  intensity  of  the  CO beam. 


VAC  " Bit 


then  the  mean  intensity  (averaged  over  the  Airy  oscillations)  as  a 


function  of  distance  is  given  by 


IL0C(X)  = ' n(x)  ±nmin 


n . = ( c/co  l) 

min  o 


In  (27),  the  factor  4 is  due  to  constructive  interference  between  the 
incident  and  reflected  CO^  fields  (this  should  really  be  considered  more 
carefully  in  the  true  cylindrical  geometry),  r)  is  the  index  of  refrac- 
tion, 


2 2 1/2 

D(x)  = (1  - u (x)/w  V' 

P o 


^min  is  t*le  nl'-n^mu,n  index  of  refraction  compatible  with  geometric  optics. 

The  maximum  possible  swelling  occurs  at  the  point  where  h(x)  = n . • 

min 

The  factor  1/2  in  equation  (27)  represents  an  average  over  the  Airy 


oscillations. 
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The  factor  ky2A2  represents  the  angular  factor,  the  factor  [ ] ^ , 

geometric  swelling,  and  the  last  factor  is  the  appropriate  frequency- 
mismatch  (or  resonant)  function  at  equal  temperatures.  If  we  intro- 
duce the  frequency-matched  wavenumber,  k , as  the  solution  to  equa- 

in 

tion  (25),  then  the  last  two  factors  in  (34)  may  be  written  as 


1 

[<3k2AD2)  + (2pkmAD)]1/2 


k 

m 

k 


2 2 
-k  /2k ^ 
m 


(35) 


U = 


c k /oj  s 1.7  (m  /m. ) 
s D o e i 


1/2 


(36) 


We  must  choose  k (the  location  of  the  frequency-matched  point  along  a 
m 

trajectory)  so  as  to  maximize  the  integrated  growth  rate.  Since  k^ky, 

and  the  exponential  must  be  prevented  from  getting  large,  we  choose 

k = k . This  allows  us  to  expand 
m y 


-k  2/2k 
m 
e 


2 


1 - ky2/2k2 


(37) 


2 2 

We  are  also  guaranteed  that  3k  / k D provided  that 


*yAD  » 2y/3 


(38) 


which  is  well  satisfied  for  the  experiment,  since  a values  are  well 
below  100.  Thus,  we  can  approximate  the  expression  in  (35)  as 


('-#)  • 


and  the  growth  rate  Y in  (34)  becomes, 

9 


where 


Yg(kx)  = IoU,p(ky4/k4)(1  " <*y2/2*2>>  . 


Xo  4cn  0 
o 


is  the  dimensionless  puiq?  intensity  at  the  point  of  maximum  Langmuir 
growth  (the  Langmuir  reflection  point) , and 


21  VAC 

o n(x  ) 
o 


(3ky2/kD2)1/2 


is  the  (swollen)  pump  intensity  in  physical  units  at  the  same  point. 

2 2 

To  recapitulate  the  meaning  of  the  factors  in  (40),  ky  /k  is  the 

A /N  2 

angular  factor  (k*y)  , k^/k  is  a pump-diminishing  factor  (due  to  geo- 
metric optics)  as  the  Langmuir  ray  passes  through  points  of  lower 

2 2 

density,  and  (k^/k) (1  - k^  /2k  ) represents  the  effects  of  imperfect 

frequency  matching,  which  become  important  away  from  the  point  k = k^. 

The  integral  for  the  integrated  growth  rate  A in  equation  (18)  may  now 

be  performed.  It  is  best  to  use  K = k A * -u>  t/2k  L as  a variable, 

x y p y 

in  place  of  t (see  equation  (14)).  Then, 
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in  order  that  K . define  the  point  at  which  I is  one-half  maximum, 
min 

Thus , 


AtK  . ) = A (0)  - 
min  z 


(47) 


It  can  be  shown  by  expansion  about  K . =0  that 

min 


A(K  . ) « A (0)  (1  - 16K  . /5TT)  , 
min  min 


(43) 


so 


K . = 

min 


ln2 


2Lk  I 

y o 


(49) 


By  invoking  equation  (15),  we  see  that 
k2 

Xmin  3(&n2)  L 1 


6x  = 3L 


4 2 ' 

(k  LI  ) 

D o 


(50) 


VI.  Now  we  can  combine  everything  to  get  a final  expression  for  a as 

a function  of  I , a,  L,  and  other  parameters  such  as  0,  r , and  <t  , which, 
o o c 

like  L and  a depend  on  the  exact  time  prior  to  peak  conpression.  The 

data  can  be  fit  quite  well  on  the  low  Iq  side  of  your  curve  (first  four 

points)  by  this  thedry,  with  L about  0.25  cm.  The  point  where 
10  2 

IQ  = 2 x 10  W/cm  requires  the  addition  of  the  nonlinear  term  in  the 
kinetic  equation.  (We  are  working  on  this  now.) 


Our  result  is  therefore 


14 


= A:  n lx  2<t>  6x6 n ( - e2A  ) , (51) 

2 2 o o c s U 

r a 
e 


where  U is  defined  by  equation  (36) - Inserting  from  (44)  and  (50) , 


with  ~ /S/a  from  (7)  , 

c 


q 

2 

r 

e 


(52) 


_1_ 

2 

a 


n 8-r  — 60 

o o a s 


3 (i.n2)  ‘ 


W! 


. « 

u 


(5TI/8)  Lk  I 
y o 


This  may  be  rewritten  as 


2 I m / 2 

r \ e / a 


(5ir/8)Lk  I 

.-V2  . Vo“V 1° 


k 2LI  2 
D o 


(53) 


2 2 2 -1 

This  is  a very  strong  function  of  temperature  since  a = kD  /k^  “ 0 , 

-3/2 

and  I “1/0  “ 0 • Thus, 

o o 


n5  CL0 
r 0 e 

_2_  <*  

2 L 


"3/2 


(54) 


As  peak  compression  is  approached,  h0  3/^2  decreases,  and  o/re  probably 
decreases,  although  the  non-exponential  factors  are  increasing. 

Now  we  put  in  parameters  which  are  independent  of  time  prior  to 
peak  compression: 


(m./m  )1/2  = 60.6  (deuterium) 

i e 

n = 9.2  x 1018  (at  the  frequency-matched  point) 
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£ = 0.25 


(spot  size) 


6fts  = di|/ ( cos  5.6°  - cos  10.4°)  = 6.53  x 10~3 

(5tr/8) Lk  I 
y o 

= 6.64  x 1017  — - . 

r 2 A^i  2 

e Do 


Parameters  associated  with  L = 0.25  cm  are 


L = 0.25  cm. 


r = 0.055  cm 
o 


(with  mask) 


0 = 0.87  keV  = 1.39  x 10~9  ergs  = 1.39  x 10_16  J 


k = 1.26  x 104  (0  = 8°) 

y s 


kD  = 1.4  x 10' 


a = k /k  = 11.1 

D y 


21.  a 

I = —ZAC — = x 8.35  x 10"14 

° 4cn  ©✓?  VAC 
o 


with  rVAC  in  units  of  W/cm  . Thus, 


a 8.67  x 103°  5,17  X 10  0lVAC 

71  = ~ e 

e VAC 


<V r (calculated) 

e 


(3.5  x 10  ) ( 13. 3)  = 4.6  x KT 


(7.2  x 1010) (311)  = 2.2  x 1013 


(5.5  x 1010) (6.4  x 102)  = 3.56  x 1013 


(1.97  x 1010)(5.2  x 104)  = 1015 


(4.4  x 1010) (1.4  x 103)  = 6 x 1013 


The  detector  mV  readings  in  the  above  chart  have  been  calculated  by 


assuming  a linear  photon  counter,  and  3 x 10  J of  scattered  energy 


for  a 110  mV  signal  strength  and  a 5J  incident  ruby  beam. 


Thus,  the  measured 


? e -?  -10 

O nd  scattered  (7.1x10  ) ( 3x10  ) 


2 2 E. 

r 4r  me 
e e 


(2. 81xl0~13) 2 


= 5.37  x 10 


for  d = 0.3  cm.  Therefore, 


V 5,371^01°  * V x 4.88  x 1011  , 


e 'measured 


where  V is  the  signal  strength  in  MV.  The  theoretical  curve  is  shown  in 


J 


Figure  4 
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APPENDIX:  Connection  of  Our  Treatment  to  Per k i ns- F 1 i ck 


I.  First,  we  show  that  the  Perkins-Flick  "threshold"  can  be  obtained 
approximately  by  our  WKB  method  under  a set  of  assumptions  appropriate 
to  their  case,  but  not  to  the  experiment: 


Assumptions : 


A.  Ignore  pun?)  swelling  and  beat  spontaneous  emission. 


B.  Assume  a well-defined  ion-acoustic  node  (requires  T » T.). 

e 1 

That  is,  Y « u>  . 

a a 

C.  Assume  the  danping  of  the  ion-acoustic  mode  is  independent  of 

k and  given  by  the  constant,  y (usually  this  is  not  the  case). 

c 

In  addition,  we  assume  a geometry  similar  to  Fig.  3,  but  with  a time- 

reversed  orbit.  Once  more,  at  the  reflection  point  k ||  E , k X Vn. 

— — o — — 

The  growth  rate  is  given  once-more  by  (22).  However,  since  we  now  have 
a well-defined  ion-acoustic  mode  (assumption  B) , the  frequency  mismatch 
function  is  given  by: 


1 1 
Im  — = Im  — 

kz  „ - — (w  -u>  +u  +iy  ) 

. C U)  L o a a 

*D2 


(Au)2  + Y 2 

cL 


where  y is  assumed  constant  (assunption  C) , and 
61 


Au)  * 0)  - U)  (k)  - c k . 

O L S 


1 


Here  c is  the  sound  speed: 
s 


c - (6  ♦ 30/^/M.1'2 

set  i 


Thus,  the  integrated  growth  rate  is  (see  equation  (22)  above) 


-/  at[Yg(kx(tl)  - irj 


E 2oj  fr  ck(k  V) 

32* "0  ~ ^ (Aoi) 2 + Y 2 ^ ' 

0 a 


where  tf  is  a final  time,  and  ky 


2 A2  = (£*e  )2.  Also,  k = (k  2(t) 


+ k 2)1//2.  With  K = k A as  a variable,  and  kx  - (u  /2L)t,  (A4) 


becomes 


, 2 ^ 
A = Ao  8 


/c  k \ y (i+k2)_i/2 

(-^L)  / dK  r 2 - Y t , (A5) 

\ Ya  / J 1+  (Aw) 2 A 

0 a 
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o = 4iTnO  ' 


. U>  t_ 
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In  order  to  do  the  integral,  we  note  that, 


..  2 2.1/2 
oj  -u)  -c  (k  +k  ) 

Aw  o L s x y 


0)  -U) 
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Intrinsic  to  the  approximation  (assumption  B)  is  that 


r i — « 1 
c k 
s y 


We  also  designate  the  frequency- matched  wavenumber  by  the  condition 

(0)  -<i>  ) » c k , so  that,  if  K = k A » then  equation  (A8)  may  be 
oLsm  m m y 

rewritten  as: 


U.kV"]  . 


so  the  integral  in  (A5)  is  of  the  form 


2.  -1/2 


/(I  + K ) ' 

~r T 

1 + r [k  - (1  + K ) 


1/2 J 2 ' 


In  this  limit  of  large  r , the  integrand  peaks  for  K ~ 1 and  K « 1 

m 

so  we  may  approximate  the  integral  is 


oo 

y . -2  K •/.  , -2  1 


1 + r — 0 1+r  — 

4 4 


/da 


The  main  contribution  to  (A5)  is  therefore 


A = A 2 /Tr  T-r-77 

o 8 r 4simr/4 
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1 

5.1 


(A14) 


If  we  define  "threshold"  as  the  value  where  A » 5,  as  do  Perkins 
and  Flick,  then  the  threshold  condition  we  find  from  (A14)  is 


o 


4Ttn0 


5.1  A(y  /c  k )** 
a s y 

Lk 

y 


A = 5 . 


(A15) 


By  comparison,  Perkins  and  Flick  find  as  their  threshold  (equation 


39,  p.  2017,  Phys.  Fluids  14  (1971): 


o 


4-rrne 


P.F. 


1.6(1  + 3T./T  ) 
l e 

k L 

y 


(A16) 


Since  T.  » T , and  V.  = 2y  (relation  between  particle  and  wave 
l e i a 

dissipation),  our  threshold  (A15)  is  higher  than  theirs  by  a factor 
of  2.25,  but  all  dimensional  factors  are  correctly  predicted.  The 
discrepancy  is  due  to  the  use  of  WKB  and,  in  particular,  arises  be- 
cause the  mismatch  function  (denominator  in  (All))  is  so  highly  peaked 
close  to  the  reflection  point  (K  = 0)  when  r » 1.  For  the  Tg  = 
case  we  treated  in  the  main-body  of  this  memo,  the  mismatch  function 
is  much  broader,  and  WKB  should  be  an  even  better  approximation'.  In 
our  treatment  we  have  the  additional  physics  of  (i)  pump  swelling 
along  the  Langmuir  ray-path,  (ii)  pump-dependent  (beat)  spontaneous 
emission  of  Langmuir  waves,  and  (iii)  ion  discreteness. 

It  is  interesting  to  note  that  if  we  take  y = c k in  (A15) , 

ci  s y 

and  v.  = 2y  , T.  = T in  (A16) , then  (A15)  predicts  a threshold 

1 d 1 G 


r 


i 
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lower  than  the  Perkins-Flick  threshold  by  a factor  of  0.56.  We  also 
note  that  if  we  take  into  account  effects  (i)  and  (iii)  above,  but 
not  beat  emission  (ii) , then  (44)  gives  a "threshold 


ELOC  > 128  ft  „ 8.1A  (Ai5) 

4iTnQ  Sir  Lk^  Lky 


which  is  within  10*  of  the  Perkins-Flick  result  (A16)  at  equal 

temperatures,  and  with  V = 2u>a.  If  spontaneous  emission  is  taken 

into  account,  then  A in  (A15)  should  be  replaced  by  A - Wn  (u/whsJ 

= A - *i£n  (oj  /c  k ),  so  the  threshold  (A15)  becomes 
o s y 


2 s.i[»-  w.  <i^Vr“>] 

r A 


LOC  __ 
4iTn0  - 


(A16) 


Lk 


y 

2A 


where  e2*  = fu  /(u  -u)  jTe2*  = I (A  measures  the  effective  amplitude 
*-  o o L 

e-folding  as  though  amplifying  from  equilibrium).  Thus,  if  A = 5, 
the  term  in  brackets  above  becomes  5 - *j&nl0  = 1.5,  and  it  is  fair 
to  say  that  beat  spontaneous  emission  lowers  the  "threshold*  by  more 
than  a factor  of  three.  Of  course,  in  terms  of  1 this  effect  is  much 
stronger,  and  without  beat  emission  we  would  have  calculated  a spec- 
trum and  cross  section  three  orders 


t 


I 


